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@ Modified DCT signal transfonnlng system. 

(57) An MDCT calculating circuit includes an xOI 
calculating circuit for multiplying input signals 
with a forward transforming window and a 
linear forward transforming unit for linear for- 
ward transforming an output signal of the cal- 
culating circuit. The linear forward 
transforming unit includes an Xo2 calculating 
circuit and an Xqs calculating circuit for pre- 
processing the output signal of the Xoi calculat- 
ing circuit and an integration and summation 
processing circuit for executing integration and 
summation processing operations on an output 
signal of the pre- processing unit. The inte- 
gration and summation processing circuit exec- 
utes an integration and summation operation 
on an N/2 number of input signals from the 
pre-processing unit by grouping a k number of 
input signals as a processing unit and iterativeiy 
executes the integration and summation pro- 
cessing operations a N/(2'^K) number of times 
for outputting a sum total of N/2 number of 
signals. 
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This invention relates to a modified MDCT signal transforming device for executing fast calculation of linear 
transformation on digital signals such as speech, audio ar picture signals. 

There has hitherto been known a transform encoding employing a so-called orthogonal transform as a 
type of the high efficiency encoding of encoding chronological sampling data such as audio signals by high- 

5 efficiency bit compression. The transform encoding means encoding after block-based orthogonal transform 
of input signals. Typical of the transform encoding is a discrete cosine transform (DCT). Such transform en- 
coding suffers from a problem of block distortion in which non-continuous block-to- block transitions are per- 
ceived as noises. The conventional practice for coping with such problem has been to have an end of a block 
overlapped with an adjacent block. According to a modified DCT (MDCT), a given block is overlapped with 

10 both adjacent blocks each by one-half a block size in such a manner as to preclude superfluous transmission 
of samples of overlapped portions. Thus the MDCT is highly suitable for high efficiency encoding. 

Such MDCT and IMDCT. which is its inverse transformation, are disclosed in "Filter Constraints in Hybrid 
MDCT Containing Plural Block Sizes" by Mochizuki. Yano and Nishlya in SHINGAKU-GIHO, CAS 90-10. 
DSP90-14, pages 55 to 60 and in "Adaptive Transform Encoding with Adaptive Block Lengths Using MDCT" 

15 by Uzu, Sugiyama, Iwatare and Nishiya in Extended Abstracts in Spring Meeting of Japan Society of Electronic 
Information Communication 1990, A-197. The above-mentioned MDCT and IMDCT are explained briefly by 
referring to Fig.1. 

Referring to Fig.1. a given block of chronological sampling data, such as Jth block, has an overlap with 
the (J-I)th block and the (J+1)th block each by a half block (50% overlap). If the number of samples of the 

20 Jth block is n, N being a natural number, there is an overlap of N/2 samples between the Jth block and the (J 
- 1)th block, while there is similarly an overlap of N/2 samples between the Jth block and the (J +1 )th block. 
Each of N samples of each of these blocks, such as the above-mentioned Jth block, is multiplied by a window 
Wh for forward transform or a pre-processing filter to produce an N number of chronological data 102. 

The characteristics of the pre-processing filter or window Wh for orthogonal transform are selected de- 

25 pending on statistic properties of input signals so that power concentration of the transformed data will be max- 
imum. By processing these N samples of the chronological data 102 by linear forward transformation by MDCT, 
an N/2 number of independent spectral data 103, that is one-half the number of input samples of the spectral 
data, are produced on the frequency domain. These N/2 number of the spectral data 103 are processed by 
inverse MDCT for producing an N number of chronological data 104 which are multiplied by a window Wf for 

30 executing inverse transform or a syntheticf liter to produce chronological data 1 05. The resulting data are sum- 
med to output results of the forward and backward blocks to restore the original chronological input data. 

Meanwhile, in practical processing operation of MDCT and IMDCT. a method for implementing fast oper- 
ation using fast Fourier transform has been proposed in A research in MDCT System and Fast Operation" by 
Iwatare. Nishiya and Sugiyama, in SHNGAK-GIHO, CAS 90-9, DSP 90-13, pages 49 to 54. The same method 

35 has also been disclosed in JP Patent KOKAI Publication No.4-44099 (1 992). The calculating method for MDCT 
and IMDCT disclosed therein is hereinafter explained. 
The definition of MDCT is given by an equation (1). 



>'o(^) =CoY>^c in) bin) COS (ic (2ic+l) (2n*N/2+l) / (2N) ) 



0 < k < N/2-1 • • ( 1 ) 

45 

In the above equation (1), x© is an MDCT input signal. N a block length, h a window function for forward 
transform, yo an MDCT output signal. Co a constant, n an integer of from 0 to N-1 . and k an integer of from 0 
to N/2-2. Since MDCT processing on sliced chronological data is executed independently for each of the blocks. 
50 the block number J is omitted. Meanwhile, since the value of Cq is not critical in carrying out the calculation 
of MDCT, the following explanation is made on the assumption that Co = 1 . For carrying out the operation. Xq 
is multiplied by a window for forward transform to find Xoi as indicated by the equation (2). 

Xoi(n) = Xo(n)h(n) 0 ^ n ^ N - 1 (2) 
From Xoi, as found from equation (2), Xo2 is calculated in accordance with equation (3). 
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-XQi(n + 3N/4) 0 < n < N/4-1 

<02(n) = I XQ^(n-N/4) N/4 < n < N-1 ••■(3) 

5 

Xo3 is calculated in accordance with equation (4). 

Xo3(n) = Xo2(2n) - Xo2(N - 1 - 2N) 0 ^ n ^ N/2 - 1 (4) 
Then, Xo3 is multiplied by a coefficient as shown by the equation (5) to produce a series of complex number 
signals Zoi. 

10 Zoi(n) = Xo3(2n)exp( - j2nn/N) 0 ^ n ^N/2 - 1 (5) 

This signal series is processed with FFT having a length of N/2 to produce a series of complex number 
signals Zo2 shown by the equation (6). 



^02^-^)=I3 {.T)exp(-j27ricn/(Ar/2) ) 

2*0 



0 < k < N/2-1 • • • (6) 



yoi is calculated from the produced series of complex number signals, as indicated by the following equa- 
tion (7). 

yoi(k) = Re(zo2(k)exp(-j27i(k + 1/2)/(2N))) 0 ^ k ^ N/2 - 1 (7) 
The above is coincident with yo defined by the equation (1 ), as demonstrated by the above-cited Research 
25 Concerning MDCT System and Fast calculation. 

On the other hand, definition of IMDCT is given by an equation (8). 



y/2-i 

x^in) =C^f{n) yAk) cos {tz {2k^l) {2n^N/2^1) / (2N)) 



0<n<N-1 •••(8) 



In the above equation (8), yi is an IMDCT input signal. N a block length, m a block number, f a window 
function for inverse transformation, x^ is an IMDCT output signal. Ci is a constant, n an integer of from 0 to N- 
1 and k an integer of from 0 to N/2-1. Since IMDCT processing on sliced chronological data is executed inde- 
pendently for each of the blocks, the block number J is omitted. Meanwhile, since the value of Co is not critical 
40 in carrying out the calculation of MDCT. the following explanation is made on the assumption that Co = 1 . For 
carrying out the calculation, y^ is re-arrayed in accordance with the following equation (9) to produce yn: 

^ yi{2k) 0 < k < N/4-1 

^ yu(^) = ' -y^(N-1-2k) N/4 < k < N/2-1 

... (9) 

50 y^i is then multiplied by a coefficient as shown by the equation (10) to produce a series of complex number 

signals z^v 

zii(k) = yii(k)exp( - j2Kk/N) O ^ k ^ N/2 - 1 (10) 
The complex number signal series is processed by inverse FFT over a length N/2 to produce a series of 
complex number signals Zi2 by the equation (11): 

55 
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.y/2-1 

^12 = X) ®^ [j2T.zik/ (N/2)} 

5 

0 < n < N/2-1 •••(11) 



The complex number signal series thus formed is multiplied by a coefficient as shown by the equation 
10 (12) and a real number part is taken out to produce x^v 

xii(n) = Re(2zi2(n) exp( - j27i(n + 1/2)/(2N)) 0 ^ n ^ N/2 - 1 (12) 
This Xii is re-arrayed, while changing its sign, as indicated by the equation (13), and is multiplied by Inverse 

n<Qfnrminn u/inHnu/ 
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20 



transforming window. 



(n) 



f {n)Xii(n+N/4) 0 S n s N/4-1 

-f (n)xll(3N/4-l-n) N/4 S n S 3N/4-1 .•.(13) 
-f (n)Xii(n-3N/4) 3N/4 S n S N/2-1 



The above is coincident with yo defined by the equation (8), as demonstrated by the above-cited Research 
Concerning MDCT System and Fast calculation. 

However, since the fast processing of MDCT for finding the spectrum for each domain having a length N 
and IMDCT which is its inverse operation are carried out with the present method using a complex number 

25 FFT having a length N/2, the N/4 log(N/2) times of multiplication of complex numbers, the N/4 log(N/2) times 
of addition of complex numbers and a work area for storing N/2 number of complex numbers are required if 
the algorithm as described above all in pp. 196 to 198 of "Fast Fourier Transform", by E. ORAN BRINGHAM, 
translated by Miyagawa and Imai, Theory of Base 2 FFT Algorithm, is used, meaning that the method cannot 
be said to be sufficiently effective if it is attempted to carry out fast processing by a small-scale processing 

30 unit. 

Various aspects of the invention are set out in claims 1 and 8 appended hereto. Preferred features are set 
out in the dependent claims. 

An embodiment of the present invention can provide a signal processing method by modified DCT in which 
the operation of the modified DCT and IMDCT, its inverse transform, may be realized by integration and sum- 
35 mation based on a K number of terms as a processing unit with the use of a small-sized processing unit. 

An embodiment of the invention provides a modified DCT signal transforming device for processing an N 
number of sample input signals, comprising a forward transforming window unit for multiplying input signals 
with a forward transforming window and a linear forward transforming unit for linear forward transforming out- 
put signals of said forward transforming window unit, characterized in that the linear forward transforming unit 
40 comprises a pre-processing section for pre-processing the output signals of the forward transforming window 
unit and an integration and summation processing section for executing integrating and summation processing 
operations on output signals of the pre-processing section, the integration and summation processing section 
executing integration and summation processing operations on an N/2 number of input signals supplied from 
the pre-processing section, the N/2 number of the input signals being previously grouped into sets each con- 
45 sisting of a K number of terms, said integration and summation operations being iteratively executed an N/(2*K) 
number of times for outputting a sum total of N/2 number of signals. 

It is noted that K in the convolution unit is a power of 2 and is not less than 4. 

The pre-processing unit outputs, with an N number of samples of input signals of from 0 to (N-1), obtained 
on multiplying the forward transforming windows in the forward transforming window section, the n + 3N/4th 
50 input signal inverted in sign as the nth intermediate signal for n of from 0 to N/4-1 . and an n - N/4th input signal 
as the nth intermediate signal for n of from N/4 to N-1. With an N number of samples of intermediate signals 
of from 0 to (N-1), the pre-processing unit outputs the 2n-th intermediate signal less the N-1-2n-th intermediate 
signal as an nth output signal. 

The integration and summation processing unit decreases the number of times of the processing opera- 
55 tions by grouping an N/2 number of input signals supplied from t he pre-processing unit into sets each consisting 
of a K number of terms, based on periodicity, for outputting the N/2 number of the signals. 

The integration and summation processing unit performs, on an N/2 number of input signals supplied from 
the pre-processing unit, a series of processing operations form of from 0 to N/(2K)-1 until the number of output 
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signals becomes equal to N/2. The series of the processing operations consists in providing, as a jth first in- 
ternnediate signal, a sum for Is of from 0 to N/2-1 of a product of the (Kl+j)th input signal multiplied by 
cos(7i(2m+1)(4K+1)/(2N). for j of from 0 to K-1 , providing, as a jth first Intermediate signal, a sum for Is of from 
0 to N/2-1 of a product of the (K(1-1)+j)th input signal multiplied by sln(7i(2m+1)(4KI+1)/(2N), for j of from K to 

5 (2K-1), providing, as a jth second intermediate signal, a sum of the jth first intermediate signal multiplied by - 
cos(2Tt(2m+1)(j/N) and the (j-K)th first Intermediate signal multiplied by -sin(27c(2m+1)(j/N), for j of from K to 
K-1. providing, as a jth second intermediate signal, a sum of the (i-K)th first intermediate signal multiplied by 
-sin(27c(2m+1)(j-K)/N) and the jth first intermediate signal multiplied by -cos(27r(2m+1)(j-K)/N), for j of from K 
to 2K-1, providing, as a k-th third intermediate signal, a sum of the jth second intermediate signal multiplied 

10 by cos(27iJk(4j+1)/(4K)), for j of from 0 to K-1, providing, as the k-th third Intermediate signal, a sum of the jth 
second intermediate signal multiplied by sin(2^k(4(j-K)+1)/(4K)), for k of from K to (2K-1), and providing, as a 
(m+kN/(2K))th output signal, a sum of the k-th third intermediate signal and a (m+kN/(2K))th third Intermediate 
signal, for k of from 0 to K-1 . 

The Integration and summation processing unit performs a processing operation of computing an output 

15 signal of the integration and summation processing unit from the second intermediate signal by a processing 
volume reduced by one half by grouping the second intermediate signals by each two terms by taking advan- 
tage of periodicity so that the number of output signals becomes equal to N/2. 

The integration and summation processing unit performs a processing operation of computing an output 
signal of the integration and summation processing unit from the second intermediate signal by a processing 

20 volume reduced by one quarter by grouping the second Intermediate signals each by four terms by taking ad- 
vantage of periodicity so that the number of output signals becomes equal to N/2. 

A modified DCT signal processing device for processing an N/2 number of independent input signals, com- 
prising a linear inverse transforming unit for linear inverse transform of input signals, and an inverse transform- 
ing window unit for multiplying an output signal of the linear inverse transform unit by an inverse transforming 

25 window, wherein the linear inverse transforming unit comprises a pre-processing unit for pre-processing an 
input signal , an Integration and summation processing unit for executing an integration and summation proc- 
essing operations on an output signal of said pre-processing unit, and a post-processing unit for post- 
processing an output signal of the integration and summation processing unit, the Integration and summation 
processing unit executing the integration and summation operations on the N/2 number of input signals sup- 

30 plied from the pre-processing unit by grouping the input signal by a K number of terms, and executing the op- 
erations Iteratively by N/(2*K) number of times for outputting a sum total of an N/2 number of signals. 
In the integration and summation processing section, K is a power of 2 which is not less than 4. 
The pre-processing unit outputs a signal equivalent to a k-th output signal which is the 2k-th input signal 
for values of k of from 0 to N/4-1 and a signal equivalent to a k-th output signal which is the (N-1-2k)th Input 

35 signal inverted in sign for values of k of from N/4 to N/2-1 . 

The integration and summation processing unit decreases the number of times of processing operations 
for integration and summation processing of N/2 numberof input signals supplied from the pre-processing unit 
based on periodicity for outputting said N/2 number of signals. 

The integration and summation processing unit performs, on an N/2 number of Input signals supplied from 

40 the pre-processing unit, a series of processing operationsfor mof from 0 to N/(2K)-1 until the numberof output 
signals becomes equal to N/2. The series of the processing operations consists in providing, as a jth first in- 
termediate signal, a sum for Is of from 0 to N/2-1 of a product of the (Kl+j)th Input signal multiplied by 
cos{7i(2m+1)(4K+1)/(2N), forj of from 0 to K-1. providing, as a jth first intermediate signal, a sum for Is of from 
Oto N/2-1 of a product of the (K(l-1)+j)th input signal multiplied by sin(7i{2m+1)(4KI+1)/(2N), forj of from Kto 

45 (2K-1). providing, as a jth second intermediate signal, a sum of the jth first intermediate signal multiplied by - 
cos(27c(2m+1)(j/N) and the (j-K)th first intermediate signal multiplied by -sin(27i(2m+1)(j/N). forj of from K to 
K-1. providing, as a jth second intermediate signal, a sum of the (j-K)th first intermediate signal multiplied by 
-sln(27r(2m+1)(j-K)/N) and the jth first intermediate signal multiplied by -cos(2;r(2m+1)(j-K)/N). for j of from K 
to 2K-1, providing, as a k-th third intermediate signal, a sum of the jth second intermediate signal multiplied 

50 by cos(27:Jk(4j+1)/(4K)), forj of from 0 to K-1, providing, as the k-th third intermediate signal, a sum of the jth 
second intermediate signal multiplied by sin(27ck(4(j-K)+1)/{4K)), for k of from K to (2K-1), and providing, as a 
(m+kN/(2K))th output signal, a sum of the k-th third intermediate signal and a (m+kN/(2K))th third Intermediate 
signal, for k of from 0 to K-1 . 

The integration and summation processing unit performs a processing operation of computing an output 

55 signal of the integration and summation processing unit from the second intermediate signal by a processing 
volume reduced by one half by grouping the second intermediate signals by each two terms by taking advan- 
tage of periodicity so that the number of output signals becomes equal to N/2. 

The integration and summation processing unit performs a processing operation of computing an output 
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signal of the integration and summation processing unit from the second intermediate signal by a processing 
volume reduced by one quarter by grouping the second intermediate signals by each four terms by taking ad- 
vantage of periodicity so that the number of output signals becomes equal to N/2. 

The pre-processing unit outputs a signal equivalent to an (n+N/4)th input signal for the value of n of from 
5 0 to N/4-1, a signal equivalent to a 3N/4-1 -nth input signal reversed in sign for values of n of from N/4 to 3N/4- 
1 and a signal equivalent to an n-3N/4th input signal reversed in sign for values of n of from 3N/4-1 to N/4-1 . 

In an embodiment of the invention, the integration and summation processing unit in the linear forward 
transforming unit or in the linear inverse transforming unit performs an integration and summation processing 
operations on an N/2 number of input signals outputted from the pre-processing unit by grouping the input sig- 
10 nals by a K number of terms and performing the operations iteratively by an N/(2*K) number of times for pro- 
ducing a sum total of N/2 signals, for simplifying the constitution. 

The invention will be described hereinafter, by way of example only, with reference to the accompanying 
drawings, in which: 

Fig.1 is a diagrammatic view showing the operating sequence of modified DCT or MDCT operations and 
15 an inverse MDCT or IMDCT operations. 

Fig .2 is a block circuit diagram showing an arrangement for realizing an MDCT signal transforming method 
according to an illustrated embodiment. 

Fig.3 is a block circuit diagram showing an arrangement for carrying out the IMDCT signal transforming 
method according to an illustrated embodiment. 
20 Fig.4 is a flow chart showing the process steps of a basic embodiment of an MDCT signal transforming 

method according to an illustrated embodiment. 

Fig.5 is a flow chart showing the process steps of a basic embodiment of an IMDCT signal transforming 
method according to an illustrated embodiment. 

Fig. 6 is a block circuit diagram showing an example of a circuit of a high efficiency encoding device for 
25 carrying out the MDCT signal transforming method according to the illustrated embodiment. 

Fig .7 is a block circuit diagram showing an example of a circuit of a high efficiency decoding device for 
carrying out the IMDCT signal transforming method according to the illustrated embodiment. 

Referring to the drawings, preferred embodiments of the signal transforming device by the modified DCT 
is explained in detail. 

30 Specifically, the signal transforming device by MDCT of the illustrated embodiment (MDCT processing cir- 

cuit 30) has an arrangement as shown in Fig.2, while the signal transforming device by IMDCT (IMDCT proc- 
essing circuit 50) has an arrangement as shown in Flg.3. Before proceeding to description of the signal trans- 
forming devices by MDCT and IMDCT of the embodiments shown in Figs.2 and 3, the basic principle of the 
signal processing method (signal transforming method) applied to these devices is explained. 

35 With the signal processing method applied to the device of the present embodiment, a coefficient table 

shown by the following equations (14) to (16) is provided for an integer m which is not less than 0 and not more 
than (N/(2K)-1): 

a^.i = cos(7i(2n + 1)(4K + 1)/(2N)) 0 ^ 1 ^ N/2 - 1 (14) 
bn,,i = sin(7i(2m + 1)(4K1 + 1)/(2N)) 
40 Cm,i = cos (27c(2m + 1 )j/(N) o ^ j ^ K - 1 (15) 

dm.i = sin(2;r(2m + 1)j/N) 
eicj = cos (27ik(4j + 1 )/(4K)) 0 ^ j ^ K/4 - 1 (16) 
fkj = sin(2tik(4j + 1)/(4K)) 
For integers m which is not less than 0 and no more than (N/(2K)-1). the following calculation is made re- 
45 peatedly. 

First, using the coefficients of the equation (14), a K numberof signals Aoi, as given by the equation (17), 
and a K number of signals Ao2t as given by the equation (18), are calculated from Xoa of the equation (4). 

^7/(2^-1 
1 = 0 



0 < j < K-1 • • • ( 17) 
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o<j<K-i -(le) 

Then, from the coefficient of the equation (15) and from these signals Aqi, A02, two signal trains or rows 
10 B01, B02. each being of a length K, are calculated as shown by the equations (19) and (20). 

Bn,.oi(i) = c^j V01G) - d^jA^,o20) 0 ^ j ^ K - 1 (19) 

Bm.02G) = - dmj Ani.oia) - C,^.J ^m^OlQ) 0 ^ j ^ K - 1 (20) 

Using the coefficients of the equation (16) and Bqi, B02, a K number of signals Xo are obtained, as shown 
by the equation (21). 

15 



X, (m^kN/ {2J0 ) e,^,B,,,, (j) f,,,B^^,, (j) 
20 ^-'^ 

0 < m < N/(2K)-1 
0<k<K-l — -(21) 

This is coincident with yo as yo defined by the equation (1), as may be seen from the following: 
First, from the relationship of the equations (1), (2), (3) and (4), the equation (22) holds: 

.v-i 



25 



30 



35 



40 



45 



50 



J2 -^01 (-^--^^Z^) cos {2tz (2k^l) (4 (n-A'/4) / (SAO ) 
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= ^ci in-N/4)cos (2- (2^:+!) (2/1-1) / (4/7) ) 



+ J2 Xo^(i3-?//4)cos(27: (2n+l) / (4iV) ) 



•v-l 

= 5^ (22-/^/4) COS (27U (2i:+l) (2n-^l) / (4/7) ) 



*V/4-l 

+ J2 ^01 (n+327/4}cos (271 (2Jc+l) (2 (ii+ic) +1) / (4N) ) 



= 53 >:oi(^-^V/4) COS {2t, {2k^l) (2n + l) / (4iV) ) 



- Xo^(n+3/7/4)cos (2?: (2JC-1) (2n-l) / (4/7) ) 



=^ -^02 COS (271 (2i:+l) (2r>+l) / (4/7) ) 



.v/2-1 

J2 ^02 (2n)cos (27: {2k^l) (4n-^l) / (4An ) 

n=0 



' X) (2n+l)cos(27i: (2i:+l) (2 (2/3+1) +1) / (4/7) ) 

x:=0 



= (2n) COS (27: (2ic+l) (4n+l) / (4/7) ) 

r.*0 
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.\V2-1 

r ^ {N'l-2n) cos (2ii (2ic^l) (2^7-1-472) / (4/7) ) 



AV2-1 

= X) -^02 (2i3)cos (2ir (2JC+1) (4^3+1) / UN) ) 



X/2-1 

- ^02 iN-l'2n) cos (27; (2;:+!) (-l-4ll) / (4N) ) 

r*0 



.V/2 -1 

= ^^02 (2n) -Xo2 {N-l-2n) ) cos (27: (2ic+l) (4n+l) / (4W) ) 

r.=0 



*V/2-l 

= ^03 (27: (2i:+l) (4n+l) / {AN) ) 

• • • (22) 

Modifying this, the equation (23) is derived. 

iV/2-l 

y^{m^kN/ {2K) ) = ^ ^oz (t: (/n+2 (kN/ {2K) ) +1) (4i-l) / (2An ) 



J2 X33(J)cos(7:(2/r?+l) (4i+l)/{2iv7) +27;ic(4i+l) MiC) 

i=0 



?r-i ii/ (2:0 -1 

= 53 X A-o3(iri+j)cos(7r (2in+l) (A(Kl^j) +1)/(2.V} 

j»0 i-O 



+ 27iic(4 {Kl^j) +1) / (4i0 ) 
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K'l /*/(2J0-l 

X) x^:, {Kl+j) COS U {2m-^l) C4A'J+1) / (2^0 

j=0 i=0 



+271 {2m+l) j/N+2Tzk(4j+l) / (4i0 ) 



=^ COS (271 (2/n+l) j7i^+27;;c(4j+l) (4/^ ) 



X iKl+j) COS (7t (2m+l) (4iCJ+l) / (2/7) ) 

i«0 



sin(27c (2;n+l) j7-V+27iic(4j + l) / (4i0 ) 



X) ^C3 (^-^+J) sin(7z (2in-l) (4JCi + l) / (2N) ) 



0<k<K-1 •••(23) 
Using the equations (14), (15), (17) and (18). the equation (24) is derived. 

{m^klJ/ {2K) ) cos (2t: (2/r?-l) j/N^27zk[ (4j +1} / {AK) Ya^^^^ (. 



sin(27: (2/n+l) j/.V+27rJc( (4j+l) / UK) )A^^^^ (j) 



=^ COS (27- (2.77+1) J /i^ COS (27X.ic(4j+l) / (4/:) ) A^,ci 

J=0 
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'Y, sin{2r.{2m^l) j/N) sin (27tic (4 j +1 ) / (4iC) )At,,oiO') 

J = 0 



K-l 

-J2 sin(27i {2in^l)j/N) cos (2u;c(4j + l) / (4in )A^^o2 U) 



cos {2n {2m^\) j / N) s±n{2^k{Aj^X) / {AK) )A^^^^[j) 



=E ^.■n,i<=os (27wic(4j+l) / (4i0 ) A^.oi (J) 



Jf-1 

■E ^;n,iSin C27t;c (4J +1) / (4iO ) A^^oi (J) 



"E ^:n.:fCOS (27u;c(4j+l) / (4jC) ) ^^^^ (j) 



■E c.t,iSin{2KJc(4j^l) /(4J0 )^^.,3,(J) 



:<-i 

= ^ cos (2TC.^c(4j^l) / i^K) ) (c,,,A,oi(J) -^..,,-.,02 (J) ) 



sin(27ric(4j>l) / UK) ) ( -d,,-A,^oi (j) "^.,^^.,02 (j) ) 



=5^cos(2irJr( (4j + l)/(4iC) )5^,oi(j) 



*E sin(2uic(4j^l) / (4i0 )B^^02O') 

J«0 



0<k<K-1 (24) 
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From the equations (19) and (20), the equation (24) is changed to the equation (25). 



0 < m <N/(2K)-1. 0 < k < K-l -..(25) 



This coincides with the equation (21) for the method embodying the present invention. It follows from this 
that the MDCT processing may be achieved with the method embodying the present invention. 

If the periodicity of the equation (16) is utilized, the amount of convolution processing of the equation (21) 
may be reduced to one half or to a quarter. 
15 That is, the coefficient of the equation (1 6) is modified to the form shown in Equations (26) and (27). 

^2k+pJ+K/2 = (2TT(2k+p) (4( j + K/2 ) + 1 )/(4K) ) 

= cos (2Tr(2k4.p) (4 j+1 )/(4K) + (2k + p)TT) 



= )^^2krpj P = 0, 1 • ■ (26) 

^2k+p,j4K/2 = (2K(2k + p) (4( j + K/2)+l )/(4K) ) 



= Sin (2TT(2k + p) (4j+i )/(4K) + (2k+p)TT) 

= ^Sk+pJ P = 0.1 - .(27) 



Therefore, as shown by equation (28), by setting, as shown by the equation (28), 

Em,o(i) = Bni.oiO) + Bm.oi(j + K/2) 
Em.ia) = B^.oia) - B^^oiO + K/2) 
35 F^,oa) = B^^oaO) + Bm.o2a + K/2) 

Fm,i(j) = B^,o2{i) - B^.o2a + K/2) 
0 ^ j ^ K/2 - 1 (28) 

Xo Is expressed as shown by the equation (29), so that the amount of integration and summation processing 
operations Is halved. 



P = 0, 1 0 < m < N/(2K)-1. 0 < k < K/2-1 •••(29) 



Besides, since the coefficient of the equation (16) may also be expressed as shown by the equations (30), 
50 (31): 



55 
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®4k+pJ+qk/4 = cos(2Tr(4k+p)(4( j + qK/4) + 1 )/(4K) ) 
= cos(2TT(4k + p)(4j+i )/(4K) + pqTr/2.): 
= cos(pqTT/2)e4^,p^j-sin(pqA/2)f4j^^p j 

p= 0,1,2,3, q= 0,1,2,3 ...(30) 

10 

^4k+p. j+qk/4 = sin(2Tr(4k+p)(4( j+qK/4)+1 )/(4K) ) 
= sin(2u(4k+p)(4j + 1 )/(4K) + pqTT/2) 
= sin (pqTr/2)e4,^,p j-sin(pqA/2)f4^,j3 j 

p = 0,1.2.3, q = 0,1,2,3 (31) 

20 the equation (32) holds. 
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j«0 Qr»0 



r/4-1 3 

to 

= E E (cos(pc77i/2) e,;,,p^^ 

15 

- s In ( pgrnc / 2 ) f^j^^^^ j ) 5^. o i O' + ^i^/ 4 ) 



r/4-i 3 

"EE (sin(pg7r/2) e^;^,^^^ 

J=0 <T*0 



25 +COS (pgTi/2) r^ic.p,^) -5^.02 



;r/4-i 3 
= J2 ^4Jt*p,:;X). (cos (pg7i/2)S,, 0,(7 +gii:/4) 

1=0 «r=0 



-sin ipqn/2) 3„^2 (j-^Q^/^) 

35 

K/4'1 3 

- E ^^k*p.jYl (sin(pg:i/2)S^^o. 



40 

+COS (pgrii/2)S^ 02 ) 
0 < m < N/(2K)-1, 0 < k < K/4-1, p = 0, 1. 2. 3 



Since one of sine and cosine values in the equation (31) is €1 witli the other value being 0, by setting so 
that, as shown in equation (33) and equation (34) 

Gm.oU) = B„.oiO) + B„,oia + K/4) + B„,oi(i + 2K/4) + M„.oi(j + 3K/4) 
G„,ia) = B„.oia) + B„.o2a + KM) - B„,oiO + 2K/4) - M„,o,G + 3K/4) 
G™,2a) = B„,oi(j) - B„.oia + K/4) - B„,oia + 2k/4) + M„,.oi(j + 3K/4) 
Gm,3a) = B„.o,G) - B„.o2a + K/4) - B„.oia + 2K/4) + M„,02(j + 3K/4) 
Gm,oO + K/4) = B„.02(j) + B„.o2(j + K/4) + B„.o2(j + 2K/4) + M„,o2(j + 3K/4) 
Gm.iO + K/4) = B™.o2a) - B,„.o20 + K/4) - B„.o2(j + 2K/4) + M„,oi(j + 3K/4) 
G,„,2a + K/4) = 8^.020) - B„,.o20 + K/4) + B™.o20 + 2K/4) - M„,o2a + 3K/4) 
G„.3(j + K/4) = B„,o2G) + B„.o,G + K/4) - B„,.o20 + 2K/4) - M„.o2(j + 3K/4) 0 S j S K/4 - 1 (33) 
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cos(2Tik(4j+1 )/(4K) ) 0 < j < K/4-1 

sin(2TTk(4( j-K/4)+l )/(4K) ) K/4 <j < K/2-1 ■••(34) 



10 



X^{in^{^k^p)N/ (210) = 9^jc.P^:}G^.p(j) 

J=0 



0 < k < K/4-1, p = 0, 1, 2, 3 •■•(35) 



t5 



holds, so that the volume of the Integration and summation operation of the equation (21) may be halved. 

Besides, In carrying out the operation of IMDCT with the signal processing method of the present inven- 
tion, coefficient tables of the equations (14), (15) and (16), which are the same as those for MDCT. are used 
for integers m which are not less than 0 and not more than (N/(2K)-1). The following operations are repeated 
20 for integers m which are not less than 0 and not more than (N/(2K)-1). 

First, using the coefficients of the equation (14), a K number of signals Poi as given by the equation (36) 

^/ (2J0-1 

25 -^'0 

0 < 0 < K-1 • - • (36) 

30 and a K number of signals P02 as given by the equation (37) 

N/ {2J0 -1 

35 i=0 

0 < j < K-1 • • . (37) 

40 are calculated from y^^ of the equation (9). 

Then, from the coefficients of equation (15) and from these signals Pqi, P02. two signal trains Q01, Q02. 
each having a length K, are calculated, as shown by the following equations (38), (39): 

Qm.OlC) = CrnjPm.Ol(j) " dmjA^.02(j) 0 ^ j ^ K - 1 (38) 

Qm,02(i) = - dr^jPrn^Olti) " C,n.An.02G) 0 ^ j ^ K - 1 (39) 

45 Then, using the coefficients of equation (16) and Bqi, B02, a K number of signals Wq are obtained, as shown 

by the equation (40): 

W,{m^kL/ (2J0 ) (j) f^,,0„_,, (j) 

J = 0 j~0 

0 < m < N/(2K)-1, 0 < k < K-1 ...(40) 

These signals Wo. obtained by the equation (40), are processed with a sign change, re-arraying and win- 
dowing for inverse transform, as shown by the equation (41), for calculating an N number of Wi: 

15 



55 



EP 0 590 790 A2 



(n)WQ(n+N/4) 0 < n < N/4-1 

W^Cn) = ( -f (n)WQ(3N/4-1-n) N/4 < n < 3N/4-1 

-f (n)WQ(n-3N/4) 3N/4 < n < N/1 

• • • (41 ) 

This may be found to be coincident with defined by equation (8) as shown by the following. 
First, an equation (42) holds from the equations (8) and (9). 

.V/2-1 

x^{n-N/A) -^C^fin) J2 ri(^)cos{2T: (Jt^l/2) (n^l/2) /N) 

Jc=0 



*V/4-l 

'-C^fin) y(2i)cos(27c(2J+l/2) (73-l/2)/AT 



.V/ 4 - 1 

- J2 y(2J+l)cos(27:(2i^l^l/2) /h^ 

1^0 



.V/4-1 

-■C^fin) y(2J)cos(27c(2i+l/2) {n^l/2) /N) 
:=r> 



.v/2-1 

- J2 y(2i+l)cos (27: (A'-1-2J) {n^l/2)/m 

1 =*V/4 
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= C^f{n) J2 y(2J)cos (27r(2i+l/2) (.^+1/2) 

- y(-^-l-2i)cos{27u(N-2i-l/2) {n^l/2) /N) 
= q/(j3} ^ y(2J)cos(27i(2J + l/2) (/3+l/2)/i^ 

*V/2-l 

"^E^ y(^-l-2i) cos (271 ("2i-l/2) (n+1/2) /iV) 
= C^f{n) J2 r(2i)cos(27i (2i+l/2) (ri + l/2)/A0 

.V/2-1 

- J2 y(-^^^-l-2i)cos (271 (2J^l/2) (n^l/2) / N) 

.V/4-1 

= qf(/7) yii(i)cos(27r (2i + l/2) (.n^l/2) /A-) 
- IT y.. (i)cos (27: (2i+l/2) {n^l/2) /N) 

*V/2-l 

-C^f{n) ^ y.. (i}cos(27r (2J-M/2) (n+1/2) /AO 

i=0 



If the equation (43) 

IJ(^) = 53 yii(^)cos(27i(2J^l/2) (i:+l/2)/^/) 

i«0 
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is set. an equation (44) 



u{n^kiV/{2K) ) =^ cos (27cJc(4j^l) /(4i0 ) O^^oiO*) 



K-l 

sin (2tuJc(4i+l) / (4iC) ) ( j) 

• • • (44) 

holds from the equations (14), (15), (36), (37), (38) and (39), in the same manner as the equation (24). 
Since the equations (45) and (46) 



uin^N) = ^ (i)cos (27C (2J-^l/2) /N) 

i«0 



= yii(J)cos(2r.(2i + l/2) (n + l/2)/iV^+(4i + l)7i) 



= - XI rii (-^) COS (27U (2i+l/2) /N) 

= -u(n) . - • (45 ) 



u(N/2-l-n)^ J2 y^^(i)cos(27c(2J+l/2) (///2 -l-n^l/2 ) 

i=0 



*V/2-l 

= X) yii(i)cos(27c(2J+l/2) {N/2-n-l/2) /N) 

1=0 



^V2-l 

= X) yii(J)cos (27U (2J+1/2) ( -Ay2 +71+1/2 ) /i^ 
= u(-N/2+n) 

= -u(N/2+n) . . . (46) 

hold from the equation (43). the equation (47) 

18 
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u{n^kN/ (2J0 ) e^,^C)„,,, (j) f-,,,0„,^^ (j) 

5 

0 < m < N/(2K)-1, 0 < k < K-1 .-.(47) 

holds on comparison with the equation (41). so that it is seen that th IMDCT calculation nnay be realized with 
10 the nnethod of the Illustrated embodiment of the present Invention. 

It is noted that, as in the case of MDCT. the volume of the integration and summation operations of the 
equation (40) may be redu ced to o ne half or to one quarter by taking advantage of t he periodicity of t he equation 
(16). 

First, if equations (48) 

^5 S^,oa) = Qm.OlO) + Qm.OlO + K/2) 

Sm.l(j) = Qm.0l(i)-Qm.0ia + K/2) 
Tm.0(j) = Qm,02(i) + Qm.oaO + K/2) 
= Qm.02(j) - Qm.02a + K/2) 

O^j^K/2-1 (48) 

20 are set, an equation (49) 

>V2-1 AV2-1 
25 J"0 

0 < m < N/(2K)-1, 0 < k < K/2-1, p = 0, 1 •••(49) 

30 holds, if the equations (26), (27) are used, so that the volume of integration and summation operations of the 
equation (40) may be reduced by one half. 

On the other hand, if the equations (30), (31) are used, an equation (50) 

35 ^^V4-l 3 

40 -sin(pqr7i/2) 0„^^ ij^gK/^) ) 

AV4-1 3 

"IT ^^k-P.jH (sin(pg7i/2)£)^^^o^(jVgi:/4) 

45 ^'^ 

+COS (pgnrc/2) ) 

0 < m < N/(2K)-1. 0 < k < K/4-1, p = 0, 1. 2, 3 

• • • (50) 

55 holds. Since one of sine and cosine values in the equation (50) is ^1, with the other value being 0, by setting 
so that, as shown in equation (51) 

Hm.oa) = Qm.oia) + Qm.oiO + K/4) + Q^,oi(j + 2K/4) + Qm^oiO + 3K/4) 
Hm.i(j) = Qm.oiO) + Qm,02a + K/4) - Q„,oiO + 2K/4) - Q^,oi(j + 3K/4) 

19 
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H^.2a) = Qm.oia) - Qm.oi(j + K/4) + Q^.otO + 2K/4) - Q^,oiO + 3K/4) 
Hm.3a) = Qm.oia) - Qm.02a + K/4) - Q^.oi(j + 2K/4) + Q^,o2(i + 3K/4) 
H„,,o(i + K/4) = Q^.o2(j) + Qm.02(j + K/4) + Q^.o2(j + 2K/4) + Q^.o2(j + 3K/4) 
Hm.iO + K/4) = Q^,o20) - Qm.oiO + K/4) - Q^.o2(j + 2K/4) + Q^.oi(j + 3K/4) 
5 H^,20 + K/4) = 0^,010) - Qm,02(j + 2K/4) + Q^.o20 + 2K/4) - Q^^oaO + 3K/4) 

Hm.3a + K/4) = Q^.o2(j) + QrT,.oi(i + K/4) - Q^,o20 + 2K/4) - Q^,oiO + 3K/4) 
O^j^K/4-1 (51) 
and by using the equation (34), an equation (52) 

10 N/7-1 

J=0 



^5 0 < k < K/4-1, p = 0, 1, 2. 3 -'(52) 

holds, so that the volume of the integration and sunnmation operations of equation (40) may be reduced to a 
quarter. 

It is seen from above that, with the method of the present invention, both MDCT defined by the equation 

20 (19 and IMDCT defined by the equation (3) may be realized by integration of a k number of terms and sum- 
mation, so that outputting of an N/2 number of signals may be executed by 

N * N/(2K) + 6N. + KN/4 
number of times of integration and summation operations and a 

K * K/2 + N + N * N/(8 * K * K) 

25 number of coefficient tables, so that the MDCT and IMDCT operations may be realized by a simpler hardware 
even although the number of times of the processing operations and the coefficient tables are increased as 
compared to the conventional method. 

Returning to Figs.2 ff., a preferred embodiment for practicing the above-mentioned principle of the illu- 
strated embodiment of the present invention is explained. 

30 That is, a signal transforming device for MDCT (MDCT calculation circuit 30) of the illustrated embodiment 

comprises an Xoi calculating circuit 33, as a forward transforming window unit for multiplying an input signal 
with a forward transforming window, and circuit components inclusive of an Xo2 calculating circuit 34 as a linear 
forward transforming unit for linear forward transforming an output signal of the xoi calculating circuit 34, and 
is adapted for processing an N number of samples of input signals. The linear forward transforming unit includes 

35 an xo2 calculating circuit 34 and an xo3 calculating circuit 35, as pre-processing units for pre-processing an out- 
put signal of thexoi calculating circuit 35, and an integration and summation circuit 40 for processing an output 
signal of the pre-processing unit with integration and summation operations. The integration and summation 
circuit 40 executes integration and summation on an N/2 number of the input signals supplied from the pre- 
processing unit. The integration is executed on each of a K number of terms of the input signals, grouped to- 

40 gether as a processing unit for integration. The integration and summation operations are iteratively performed 
an N/(2*K) number of times for outputting a sum total of N/2 signals. 

In a flow chart of Fig.4, the signal transforming sequence in the embodiment of Fig.2 for executing MDCT 
operations by the method of the present invention is shown. At a first step SI in Fig.4, chronological sample 
data, such as PCM audio data, are grouped into blocks each made up of a predetermined number of. herein 

45 N samples. The blocks are set so that an overlap with neighboring blocks amounts to 50%, that is the neigh- 
boring blocks are overlapped by N/2 samples, as described in connection with Fig.1, and sample data of the 
respective blocks are multiplied by a forward transforming window Wh, as described in connection with Fig.1. 
These filtered N-samples of the chronological data Xoi are modified at the next step S2 as shown by the equa- 
tions (3) and (4) for transformation into an N/2 number of real number data X03. At the next step S3, calculations 

50 shown by the equations (17) to (20), (29) and (31) are executed on X03 to find the N/2 number of spectral data 
Xo which are outputted as MDCT output data. 

The above-described transformation for MDCT is executed by a signal transforming device for MDCT of 
the present embodiment, an illustrative hardware for realizing the transformation for MDCT, shown in Fig.2. In 
the signal transforming device, shown in Fig.2, chronological sample data Xq, stored in chronological sample 

55 buffer 32, is transformed by an MDCT calculating circuit 30 for executing the signal transforming method for 
MDCT into spectral data which is outputted at an output terminal 39. 

Specifically, the MDCT calculating circuit 30 comprises an Xoi calculating circuit 33 for executing a proc- 
essing on the chronological sample data Xq according to equation (2) to find the chronological data Xqi, an X02 
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calculating circuit for executing calculations according to equation (3) on the chronological data xoi to find data 
Xo2, an calculating circuit 35 for executing calculations according to equation (4) on the data Xo2 to find data 
Xo3, and an integration and summation circuit 40 for executing calculations according to equations (17) to (20), 
equations (33) and (35) on the data X03 to find spectral data Xq. 

5 Specifically, the integration and sumnrjation circuit 40 connprisesan Am,oi calculating clrcuit41 for executing 

a processing according to equation (1 7) using data X03 and an output of an X calculating circuit 46 as later ex- 
plained to find data A^.oi* an A^.oa calculating circuit 42 for executing a processing according to equation (18) 
using the data X03 and the output of the X calculating circuit to find data Am,o2. a Bm,oi calculating circuit 43 for 
executing a processing according to equation (19) to find data 8^,01. a calculating circuit 44 for executing 

10 a processing according to equation (20) to find data Bm,o2> a calculating circuit 45 for executing a processing 
according to equation (33) using the data 8,^,01 and data 8^,02 to find data G^, and an Xq calculating circuit 46 
for executing a processing according to equation (35) to find spectral data Xq. 

Meanwhile, the Xoi calculating circuit 33 as a forward transforming window unit, and the X02 calculating 
circuit 34, as a pre-processing unit for the linear forward transforming unit, may be grouped into one circuit. 

15 A signal transforming device for IMDCT according to an embodiment of the present invention (IMDCT cal- 

culating circuit 50) comprises a linear inverse transforming unit for executing linear inverse transform on an 
input sig nal and an inverse transforming window unit for multiplying an output signal of the linear inverse trans- 
forming unit with an inverse transforming window, as shown in Fig.3, and is adapted for performing inverse 
transform operations on an N/2 number of independent input signals. The nonlinear inverse transforming unit 

20 comprises a y^ calculating circuit 52 for pre-processing an input signal, an integration and summation circuit 
60 for executing integration and summation on an output signal of the y^^ calculating circuit 52, and a Wi cal- 
culating circuit 56 for post- processing an output signal of the integration and summation circuit 60. The inte- 
gration and summation circuit 60 executes integration and summation operations on the N/2 number of input 
signals supplied from the y^^ calculating circuit 52 by grouping the input signal into groups each made up of K 

25 terms, and iteratively performs the operations N/(2*K) times to output a sum total o N/2 number of signals. 
Meanwhile, in the example of Fig. 3, the Wi calculating circuit 56 is composed of the inverse transforming unit 
and a post-processing unit of the linear inverse transforming unit. 

In a flow chart of Fig.5, the signal transforming sequence in the embodiment of Fig.3 for executing IMDCT 
operations by the method of the present invention is shown. At a first step S11 shown in Fig. 5, an N/2 number 

30 of spectral data y^ is processed by sign exchange and re-arraying as shown by equation (9) for being converted 
into an N/2 number of real number data y^^ At the next step SI 2, the real number data yi^ are processed by 
processing operations in accordance with equations (32) to (35), (43) and (44) to generate an N/2 number of 
data Wq. At the next step S15, the data Wq are processed by transform processing as shown by equation (37) 
to find an N number of real number data Wi . 

35 The above-described transform operations for IMDCT and the transform processing for constructing 

chronological data from the output of the IMDCT transform operations are executed by a signal transforming 
device forMDCTofthe present embodiment, an illustrative hardware for realizing the transformation for IMDCT, 
shown in Fig.3. The spectral data yi from input terminal 51 are multiplied by an inverse transform window by 
a chronological data buffer 57 and subsequently summed to IMDCT output data of both neighboring blocks by 

40 an overlap addition circuit 58 before being outputted as chronological sample data at output terminal 59. 

Specifically, the IMDCT calculating circuit 50 comprises a y^ calculating circuit 52 for executing calcula- 
tions on the spectral data yi in accordance with the equation (9) to find data y^, an integration and summation 
circuit 60 for executing calculations on data y^i in accordance with the equations (36) to (39), (51) and (52) to 
find data Wq. and a Wi calculating circuit 56 for executing calculations on data Wq in accordance with the equa- 

45 tion (41) to find data 

On the other hand, the integration and summation circuit 60 comprises a Pm,oi calculating circuit 61 for 
executing calculations on the data y^^ and an output of a Wq calculating circuit 66 as later described to find 
data Pm,oi. a Pm,o2 calculating circuit 61 for executing calculations on the data y^ and the output of Wq calcu- 
lating circuit 66 as later described to find data Pm,o2. a Qm,oi calculating circuit 64 for executing calculations on 

50 the data Pm,oi and Pm,o2 in accordance with the equation (38) to find data Qm.oi. a Qm,o2 calculating circuit 64 
for executing calculations on the data Pm^oi and Pm,o2 in accordance with the equation (39) to find data Qm,02. 
an Hm calculating circuit 65 for executing calculations on data Qm,oi and Qm,o2 in accordance with equation (51) 
to find data Hm and a Wq calculating circuit 66 for executing calculations on the data in accordance with 
equation (52) to find data Wq. 

55 Referring to Fig. 6. explanation is made of a concrete example of a high efficiency encoding device, to which 

MDCT as one of the above-described signal transforming operations of the illustrated embodiments is applied. 
The concrete high efficiency encoding device shown in Fig. 6 makes use of techniques of sub-band coding, 
adaptive transform coding and adaptive bit allocation, That is, with the high efficiency encoding device, shown 
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in Fig.6, input digital signals, such as input PCM audio signals, are split into plural frequency bands so that the 
bandwidth becomes broader towards the higher frequency. The band-limited signals are processed by MDCT, 
a kind of orthogonal transform processing operations, for each frequency band, to generate frequency-domain 
spectral data, which are then encoded by adaptive bit allocation for each of so-called critical bands. 

5 That Is. referring to Fig.6, audio PCM signals in e.g. a range of from 0 to 20 kHz are supplied to input ter- 

minal 11. The input signals are split by a band-splitting filter 12. such as a so-called QMF filter, into a band of 
0 to 10 kHz ad a band of 10 to 20 kHz. the band of 0 to 10 kHz is split by a similar band-splitting filter 13, such 
as QMF filter, into a band of 0 to 5 kHz and a band of 5 to 10 kHz. The signals in the band of from 10 kHz to 
20 kHz from band-splitting filter 13 are transmitted to the modified DCT (MDCT) circuit 14 of the illustrated 

10 embodlmentfor MDCT processing operations. Similarly, the signals in the band of 5 to 10 kHz and in the band 
of 0 to 5 kHz, supplied from the band-splitting filter 1 3, are transmitted to MDCT circuits 15. 16, respectively, 
for MDCT processing operations. 

The spectral data on the frequency domain or coefficient data, generated by MDCT processing operations 
by the MDCT circuits 14 to 16, are grouped for each of the critical bands so as to be supplied to an adaptive 

15 bit allocation and encoding circuit 17. 

Meanwhile, the critical band means a frequency band splitting system which takes human acoustic sense 
into consideration. That Is, the critical band Is a band of a narrow bandwidth noise which masks a pure tone 
having a frequency in the vicinity of that of the noise. The critical band becomes broader in bandwidth towards 
a higher frequency. The entire audible frequency range of from 0 to 20 kHz is split Into e.g. 25 critical bands. 

20 The adaptive bit allocation and encoding circuit 1 7 normalizes each of the spectral signals by a scale factor, 

that is a maximum value of the absolute values of the spectral signals contained In each critical band, from 
one critical band to another. The circuit 17 then re-quantlzes the normalized spectral signals with a sufficient 
number of bits to mask the quantization noise by the critical band signals to output the re-quantized spectral 
signals along with the scale factor as found for each critical band and the number of bits employed for requan- 

25 tizatlon. The data encoded in this manner are outputted via output terminal 18. 

Referring to Fig.7. explanation is made of a concrete example of a high efficiency decoding device, to which 
IMDCT as one of the above-described signal transforming operations of the illustrated embodiments Is applied. 

Referring to Fig.7. the above-mentioned scale factor, the number of bits employed for requantizatlon and 
requantized and encoded spectral signals are entered at input terminal 21 and spectral signals are generated 

30 by a spectral signal decoding circuit 22 from these Input data. Of these spectral signals, the spectral signal In 
the ranges of 10 to 20 kHz, 5 to 10 kHz and 0 to 10 kHz are transformed into signal waveform data of the 
respective bands by IMDCT circuits 23. 24 and 25, respectively. Of these signal waveform data in the three 
frequency bands, the signal waveform data in the band of 0 to 5 kHz are combined with the signal waveform 
data in the band of 5 to 10 kHz by a band combining circuit 26 into signal waveform data In the band of 0 to 

35 10 kHz which is then combined with the signal waveform data in the band of 20 to 20 kHz by a band combining 
circuit 27 into signal waveform data over the entire frequency range which Is outputted at output terminal 28. 

It should be noted that the present invention is not limited to the above-described embodiments. For ex- 
ample, the device of the present invention may be applied to a variety of transform encoding circuit or a decoder 
In addition to the high efficiency encoding/decoding device shown in Flgs.6 and 7. 

40 The forward transform calculations for modified DCT (MDCT) may be performed efficiently on an n number 

of chronological sample data sliced with 50% overlap with both neighboring blocks, using FFT with a length 
of N/4, for reducing the processing volume or work area for constructing a high efficiency encoding device for 
audio signals. Besides. Inverse transform calculations for modified DCT (MDCT) may be performed efficiently 
on an N/2 number of independent spectral data using FFT with a length of N/4 for reducing the processing 

45 volume or work area for constructing a high efficiency decoding device for audio signals. 



Claims 

50 1. A modified DCT signal transforming device for processing an N number of sample Input signals, compris- 
ing a forward transforming window unit for multiplying input signals with a forward transforming window 
and a linear forward transforming unit for linear forward transforming output signals of said forward trans- 
forming window unit, characterized In that said linear forward transforming unit comprises a pre-process- 
ing section for pre-processing the output signals of said forward transforming window unit and an inte- 

55 gratlon and summation processing section for executing Integrating and summation processing operations 

on output signals of said pre-processing section, said integration and summation processing section exe- 
cuting integration and summation processing operations on an N/2 number of input signals supplied from 
said pre-processing section, said N/2 number of the input signals being previously grouped Into sets each 
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consisting of a K number of terms, said integration and summation operations being iteratively executed 
an N/(2*K) number of times for outputting a sum total of N/2 number of signals. 

The modified DCT signal transforming device as claimed in claim 1 wherein Kin said integration and sum- 
mation processing section is a power of 2 which is not less than 4. 

The modified DCT signal transforming device as claimed in claim 1 wherein, with an N number of samples 
of input signals of from 0 to (N-1 ), obtained on multiplying the forward transforming windows in the forward 
transforming window section, the pre-processing unit outputs the n+3N/4th input signal inverted in sign 
as the nth intermediate signal for n of from 0 to N/4-1 , and an n-N/4th input signal as the nth intermediate 
signal for n of from N/4 to N-1, and wherein, with an N number of samples of intermediate signals of from 
0 to (N-1), the pre-processing unit outputs the 2n-th intermediate signal less the N-1-2n-th intermediate 
signal as an nth output signal. 

The modified DCT signal transforming device as claimed in claim 1 wherein the integration and summation 
processing unit decreases the number of times of the processing operations by grouping an N/2 number 
of input signals supplied from said pre-processing unit into sets each consisting of a K number of terms, 
based on periodicity, for outputting said N/2 number of the signals. 

The modified DCT signal transforming device as claimed in claim 1 wherein said integration and summa- 
tion processing unit performs, on an N/2 number of input signals supplied from said pre-processing unit, 
a series of processing operations for m of from 0 to N/(2K)-1 until the number of output signals becomes 
equal to N/2, said series of the processing operations consisting in providing, as a jth first intermediate 
signal, a sum for Is of from 0 to N/2-1 of a product of the (Kl+j)th input signal multiplied by 
cos(n(2m+1)(4K+1)/(2N), for j of from 0 to k-1, providing, as a jth first intermediate signal, a sum for Is of 
from 0 to N/2-1 of a product of the (K(l-1)+j)th input signal multiplied by sin(n(2m+1)(4KI+1)/(2N), for j of 
from K to (2K-1), providing, as a jth second intermediate signal, a sum of the jth first intermediate signal 
multiplied by -cos(2n(2m+1 )(j/N) and the 0-K)th first intermediate signal multiplied by - sin{2;r(2m+1 )0*/N), 
for j of from K to K-1 , providing, as a jth second intermediate signal, a sum of the (j-K)th first intermediate 
signal multiplied by -sin(27c(2m+1){j-K)/N) and the jth first intermediate signal multiplied by -cos(2?t(2m+1)(j- 
K)/N), for j of from K to 2K-1, providing, as a k-th third intermediate signal, a sum of the jth second inter- 
mediate signal multiplied by cos(27iJk(4j+1)/(4K)), for j of from 0 to K-1 , providing, as the k-th third inter- 
mediate signal, a sum of the jth second Intermediate signal multiplied by sin(27ik(4(j-K)+1)/(4K)), for k of 
from K to (2K-1), and providing, as a (m+kN/(2K))th output signal, a sum of the k-th third intermediate 
signal and a (m+Kn/(2K))th third intermediate signal, for k of from 0 to K-1. 

The modified DCT signal transforming device as claimed in claim 5 wherein said integration and summa- 
tion processing unit performs a processing operation of computing an output signal of the integration and 
summation processing unit from the second intermediate signal by a processing volume reduced by one 
half by grouping the second intermediate signals by each two terms by taking advantage of periodicity 
so that the number of output signals becomes equal to N/2. 

The modified DCT signal transforming device as claimed in claim 5 wherein said, integration and sum- 
mation processing unit performs a processing operation of computing an output signal of the integration 
and summation processing unit from the second intermediate signal by a processing volume reduced by 
one quarter by grouping the second intermediate signals by each four terms by taking advantage of per- 
iodicity so that the number of output signals becomes equal to N/2. 

A modified DCT signal processing device for processing an N/2 number of independent in put signals, com- 
prising a linear inverse transforming unitfor linear inverse transform of input signals, and an inverse trans- 
forming window unit for multiplying an output signal of said linear inverse transform unit by an inverse 
transforming window, 

said linear inverse transforming unit comprising a pre-processing unitfor pre-processing an input 
signal, an integration and summation processing unit for executing an integration and summation proc- 
essing operations on an output signal of said pre-processing unit, and a post-processing unit for post- 
processing an output signal of said integration and summation processing unit, said integration and sum- 
mation processing unit executing the integration and summation operations on the N/2 number of input 
signals supplied from said pre-processing unit by grouping the input signal by a K number of terms, and 
executing the operations iteratively by N/(2*K) numberof times for outputting a sum total of an N/2 number 
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of signals. 

The modified DCT signal transforming device as claimed in claim 8 wherein K in said integration and sum- 
mation processing section is a power of 2 which is not less than 4. 

The modified DCT signal transforming device as claimed In claim 8 wherein said pre-processing unit out- 
puts a signal equivalent to a k-th output signal which is the 2k-th input signal for values of k of from 0 to 
N/4-1 and a signal equivalent to a k-th output signal which is the (N-1-2k)th input signal inverted in sign 
for values of k of from N/4 to N/2-1 . 

The modified DCT signal transforming device as claimed in claim 8 wherein said integration and summa- 
tion processing unit decreases the number of times of processing operations for integration and summa- 
tion processing of N/2 number of input signals supplied from said pre-processing unit based on periodicity 
for outputting said N/2 number of signals. 

The modified DCT signal transforming device as claimed in claim 8 wherein said integration and summa- 
tion processing unit performs, on an N/2 number of input signals supplied from said pre-processing unit, 
a series of processing operations for m of from 0 to N/(2K)-1 until the number of output signals becomes 
equal to N/2, said series of the processing operations consisting in providing, as a jth first intermediate 
signal, a sum for Is of from 0 to N/2-1 of a product of the (Kl+j)th input signal multiplied by 
cos(7c(2m+1)(4K+1)/(2N), for j of from 0 to K-1, providing, as a jth first intermediate signal, a sum for Is of 
from 0 to N/2-1 of a product of the (K(l-1)+j)th input signal multiplied by sin(7i(2m+1)(4KI+1)/(2N), for j of 
from K to (2K-1), providing, as a jth second intermediate signal, a sum of the jth first intermediate signal 
multiplied by -cos(27r(2m+1)(j/N) and the (j-K)th first intermediate signal multiplied by - sin(2n(2m+1 )0/N), 
for j of from K to K-1, providing, as a jth second intermediate signal, a sum of the (j-K)th first intermediate 
signal multiplied by -sin(27i(2m+1 )0-K)/N) and the jth first intermediate signal multiplied by -cos(27i(2m+1)(j- 
K)/N), for j of from K to 2K-1, providing, as a k-th third intermediate signal, a sum of the jth second inter- 
mediate signal multiplied by cos(2iiJk(4j+1)/(4K)), for j of from 0 to K-1. providing, as the k-th third inter- 
mediate signal, a sum of the jth second intermediate signal multiplied by sin(27ik(4{j-K)+1)/(4K)), for k of 
from K to (2K-1), and providing, as a (m+kN/(2K))th output signal, a sum of the k-th third intermediate 
signal and a (m+kN/(2K))th third intermediate signal, for k of from 0 to K-1. 

The modified DCT signal transforming device as claimed in claim 12 wherein said integration and sum- 
mation processing unit performs a processing operation of computing an output signal of the integration 
and summation processing unit from the second intermediate signal by a processing volume reduced by 
one half by grouping the second intermediate signals by each two terms by taking advantage of periodicity 
so that the number of output signals becomes equal to N/2. 

The modified DCT signal transforming device as claimed in claim 12 wherein said integration and sum- 
mation processing unit performs a processing operation of computing an output signal of the integration 
and summation processing unit from the second intermediate signal by a processing volume reduced by 
one quarter by grouping the second intermediate signals by each four terms by taking advantage of per- 
iodicity so that the number of output signals becomes equal to N/2. 

The modified DCT signal transforming device as claimed in claim 8 wherein said pre-processing unit out- 
puts a signal equivalent to an (n+N/4)th input signal for the value of n of from 0 to N/4-1 , a signal equivalent 
to a 3N/4-1-nth input signal reversed in sign for values of n of from N/4 to 3N/4-1 and a signal equivalent 
to an n-3N/4th input signal reversed in sign for values of n of from 3N/4-1 to N/4-1 . 
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